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Abstract 

An attempt has been made to combine a wave solu- 
tion method and an unsteady flow computation to pro- 
duce an integrated aeroacoustic code to predict far-field 
jet noise. An axisymmetric subsonic jet is considered 
for this purpose. A fourth order space accurate Pade 
compact scheme is used for the unsteady Navier-Stokes 
solution. A Kirchhoff surface integral for the wave equa- 
tion is employed through the use of an imaginary sur- 
face which is a circular cylinder enclosing the jet at a 
distance. Information such as pressure and its time and 
normal derivatives is provided on the surface. The sound 
prediction is performed side by side with the jet flow 
computation. Retarded time is also taken into consid- 
eration since the cylinder body is not acoustically com- 
pact. The far-field sound pressure has the directivity 
and spectra show that low frequency peaks shift to- 
ward higher frequency region as the observation angle 
increases from the jet flow axis. 

Introduction 

Jet noise suppression has appeared as a critical is- 
sue for the viability of future supersonic flight. The FAR 
36 Stage III imposes the same noise limitation on future 
supersonic commercial flight as it does on subsonic air- 
craft. This brings forth a challenging task to reduce the 
jet noise in the High Speed Civil Transport Program. 
In order to accomplish this task, assessment of far-field 
noise generated by a jet plume is a prerequisite to the 
design of the engine. Jet noise is generated as a byprod- 
uct of the plume flow behind the exhaust nozzle. Flow 
turbulence has been believed to be a source of the sound. 
Lighthill 1 showed that the flow turbulence, which is re- 
ferred to as the fluctuating Reynolds stress or LighthilPs 
tensor, is the sound source. Since then, turbulent flow, 
which had been generally accepted as a totally chaotic 
entity, has been a central theme in the aeroacoustic re- 
search. For an estimation of the sound pressure using 
the acoustic analogy approach, the two-point fourth or- 
der correlation of the fluctuating Reynolds stress must 
be computed. To make this tractable, Proudman pur- 
sued a noise generation theory that assumes the isen- 
tropic turbulence. A variety of manipulation and mod- 
elling of the flow turbulence has been made based on a 
fully turbulent assumption since then. 

Freymuth 3 observed organized large eddy struc- 
tures in a separated flow of a jet. Brown and Roshko 4 


also found large vortical structures in a free shear layer. 
These findings of the vortical pattern in the free shear 
layer filled the gap between an initial wave region, in 
which a linear theory is applied, and the fully turbulent 
downstream region. The flow regime dominated by the 
large vortical structure is not fully random and is pre- 
dictable in a deterministic way. This organized structure 
maintains its identity up to the point where the poten- 
tial core begins to collapse but is still discernable even in 
the fully turbulent region far downstream. Winant and 
Browand 5 reported that a mechanism of the mixing layer 
growth is an interaction of adjacent large vortices. These 
investigators have shown that the flow in free shear lay- 
ers such as jet and plane mixing flow is well behaved and 
more organized than previously thought. Shear flow is 
dominated by large vortical structures, which are very 
predictable and controllable. This shear flow, which had 
been thought to be fully turbulent and therefore random 
and chaotic, has become research subject with a quite 
different perspective since the observation of these or- 
ganized structures. A decomposition of the fluctuating 
flow quantity into the organized flow entity and the fully 
random entity makes it possible to study turbulent shear 
flows in a certain deterministic way. 

Experiments 6 have shown that the sound power 
emitted from the jet column is greatest within 4 or 5 di- 
ameters downstream, and then decays rapidly through a 
transition region. This indicates that the initial develop- 
ment of the jet, before it becomes fully turbulent, should 
be clearly resolved so that an accurate noise prediction 
can be made. This region is characterized by large vor- 
tical structures and is not fully turbulent, which gives 
the motivation that we solve the unsteady flow equation 
directly to provide the sound source for an acoustic com- 
putation of the far-field noise. Numerical solution of tur- 
bulent flow is difficult because the turbulent flow field is 
made up of a range of length scales from the Kolmogorov 
scale to the integral scale. If numerical mesh size can 
be made fine enough to resolve the smallest scales which 
dissipate the kinetic energy, then direct numerical sim- 
ulation (DNS) is the tool to obtain the entire turbulent 
flow structure. However, the dissipative scale becomes 
finer as the Reynolds number is increased and practical 
hardware limitations are rapidly reached. Therefore, the 
DNS method is limited to simulating only low Reynolds 
number turbulence. For practical computation of higher 
Reynolds number flows, small scale fluctuations can be 
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modeled so that desired large scale eddies can be com- 
puted directly, while proper dissipation is provided by 
the small scale eddy model. This approach, which is 
referred to as large eddy simulation (LES), has been 
successfully employed in many flows with practical ap- 
plications. 

In order to obtain the flow field as the source of 
sound using DNS or LES, the simulations must be per- 
formed using numerical techniques with minimal distor- 
tion and diffusive characteristics. The source of numer- 
ical diffusion and phase error is known to be mainly 
from the numerical formulation of the convective terms. 
These numerical artifacts get worse for high Reynolds 
number flow simulations. Typically, free shear flows of 
interest have very high Reynolds numbers. Therefore, a 
higher order accurate numerical scheme which meets the 
previously mentioned requirements is needed. Fourth 
order Pade compact differencing scheme with a disper- 
sion relation preserving property is used here. 

It is the purpose of this paper to present a method 
to predict the far-field pressure directly from the numeri- 
cally generated unsteady flow solution without recourse 
to empirical factors. Therefore, only an axisymmetric 
laminar jet is considered in the process of incorporat- 
ing a wave solution into the higher order accurate flow 
solver. Furthermore, the flow solution is limited to the 
subsonic case since supersonic jets often generate shock 
related noise in addition to the shear noise due to flow 
turbulence, which makes the problem more complicated. 
The Kirchhoff surface integral method is chosen for the 
solution of the wave equation. The acoustic result ob- 
tained is the far-field sound caused by the wavy motion 
and large vortical structure of the free shear layer. Fine 
scale random turbulence is not addressed in this study. 

Governing Equation of Fluid Flow 
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where superscript T is the transpose of the matrix, (u, v) 
are the velocities in (i,r), and the source term s and 
stress tensors are defined as: 

S = - - [g + (0,0,p- r^,0) T ] 
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du dv 

r xx = AV ■ u + 2p— , r rr = AV • u + 2p— , 


dx’ 


dudv 

Txr ~ M ^dr + dx^ : ^ = AV u + 2 ^ r 


dr 

V 


Heat conductivity k and viscosity p are scaled by k* and 
p * , which are the values at T* . The second viscosity A 
is — | p. Dimensionless numbers, Re and Pr, are the 
Reynolds and Prandtl numbers defined to be p*u*l*/p * 
and p*c p /kt, respectively. The dimensionless speed of 
sound c becomes y/T/M r . Equation (1) is written in 
generalized coordinates £ and 77 as: 


The variables T, p, x, u,*,p, and e are dimensionless 
quantities of temperature, density, position, velocity, 
time, pressure, and total energy per unit mass based on 
reference quantities T* , p * , /* , u* , t* , p* and e* , respec- 
tively. Additional definitions are t* — l r /u T) p * = p T u 
and e* — v%. Then the equation of state becomes the 
following: 

P = ^2 with T = T(7 - 1 )M r 2 (e - y - ) 

where 7 is the ratio of specific heats, M r = u*/y/jRT * , 
here R is the gas constant. The Navier-Stokes equa- 
tions for axisymmetric flow in the cylindrical coordinates 
(x,r) are written as : 
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where Q = J q, 


S = Js (2) 


J = (ZsVr-Zr n X )~\ F = J^f-Krg), G = J (r, x i+T, r g) 


Formulation of Difference Scheme 

To obtain a first derivative of /(x), a Pade compact 
differencing 7,8,9 is formulated as : 
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where h is the mesh size. If we use the truncated Taylor 
series to make equation (3) a fourth order approxima- 
tion, the values of a and b are : 
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Equation (3) can render a sixth order accuracy if a is 
set to |. In a general fourth order formulation, a is a 
free parameter and it will be optimized by the disper- 
sion relation preserving concept proposed by Tam and 
Webb 10 as described below. Fourier transform of f(x) 
and its inverse are defined as : 


accurate and the intermediate variables are not stored 
at every stage. 

A numerical dissipation term D is added during 
the fourth stage to enhance the numerical stability. The 
dissipation term is introduced to be of sixth order so 
that our fourth order accuracy remains intact. 


>^wJZ h ^"’ dw 


< d 6 q 3 6 q>. 
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where w e is a constant and <9 6 q/d£ 6 is given by : 


Fourier transforming equation (3) by the above defini- 
tion gives : 

i (1 + 2acos«)/(u/) = i (asinK + ^sin 2k)}{w) 

This expression indicates that the wave number is de- 
formed by our discretization process into a different 
wave number as : 


£6 

— = 15 (qi + ij + q*-ij) - 6(q,-+2j + Qi-2;) 

+ (q*+3j 4- q«-3j) — 20q : ; 


The derivative 3 6 q /drf in the 77 direction is obtained 
in the similar manner. This numerical dissipation is 
applied to internal points. 

Boundary Condition 


a sin k 4- |sin 2 k 

k = - t 0 2 ( 4 ) 

1 -f 2ct COS K 

where k is the input wave number defined to be wh and 
k the deformed response wave number. In reference [10] 
the optimum ot is chosen so that the following function 
K is minimized, (i.e dK/dot = 0) 

r 17 

K = I (/c — k) 2 dn 

J-tc/2 

The optimum value of ot is computed to be 0.35619. Fig- 
ure 1 shows the wave relation. The straight line is for the 
exact derivative. It is also found that the fourth order 
approximation with the optimum value of a has better 
wave performance than the sixth order approximation 
in the most compact form. 

The four-stage Runge-Kutta technique 11 is adopted 
for an explicit time advancement formulation. To obtain 
new flow variables at t = (n 4- 1) At from known data at 
t = nAt, equation (2) is used to advance the solution in 
time as follows : 

QO) _ Q" = ai AtW< 0) 

Q( 2 )-Q n = a 2 A tW* 1 * 

Q( 3 ) - Q n = a 3 AiW (2 > 
q( 4 ) _ Q n = a 4 AtW^ + D 

where denotes S - dT /d£ - dG/drj evaluated at 

the jb-th stage. The stage 0 and 4 are at the time, nAt 
and (n+l)At . The parameters, ai ,02,^3, c* 4 , are given 
to be This time difference is second order 


The boundary treatment considered here is a com- 
bination of characteristic and algebraic boundary condi- 
tions. The characteristic boundary condition solves the 
governing equation in a characteristic form in each coor- 
dinate direction and the algebraic boundary conditions 
are the given boundary conditions such as temperature, 
total temperature, velocity, etc. Equation (2) is written 
in a non-conservative form as: 


. 3q 5q 

m +A di + B ^- Sr 


where s r = r _1 (q x F - f r G + s), A = 5Fi/i9q,B — 
dGi/dq , and the subscript i denotes the inviscid part. 
For the £ direction, equation (2) with a transformation 
dq = R^dq becomes 

f + R r 'AR ( | + R t -‘BR ( fi = 0 

If we construct the matrix such that the eigenvectors 
of the matrix A constitute its columns then the matrix 
by a similarity transform becomes a diagonal matrix, 
whose entries are the eigenvalues of A such that 


R^ 1 ARf = A { = diag (U, U,U + a^,U - a ? ) 


where a e = Cyjg+tf. Here, c is the speed of sound 

defined to be yj T/M r , and U — £ x u + £yV- In fhe same 
way, for the tj direction the diagonal matrix becomes . 

R^BRq = A n = diag(V, V, V + a„, V - a„) 
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where a n = CyJrjl + f7 2 , and K = 7? r t£-f r) y v. The charac- 
teristic equations are then rewritten in each coordinate 
direction to be: 
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1^5 , * p -i. w ^ | 

* m + A? Rf ** + 
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Nonreflecting boundary conditions 12,13 can be con- 
structed by setting any eigenvalue, which is the element 
of A, to be zero, if a wave is incoming towards the com- 
putational domain. 


Jet Flow Calculation 


The subsonic jet Mach number is 0.6 and the ambi- 
ent air is at M=0.2 and the two streams are brought to 
be mixed at the same temperature. The reference length 
l* is taken to be the nozzle radius R. u * is the average 
velocity of the two streams and M r is the average Mach 
number of the two streams. Since the temperature vari- 
ation is assumed small over the entire flow domain the 
viscosity and the heat conductivity are held constant 
at the reference temperature T* = 298° K. The inlet 
boundary conditions are given by : 


u — \ — A s tanh(20(r - 1)), v — 0, T = 1 (7) 


and there is one characteristic equation for the outgo- 
ing wave. The shear ratio A, is defined to be A u* /u*, 
where A u* is the velocity difference in the two streams. 
The characteristic boundary conditions at the exit and 
side boundary planes are given as described in equations 
(5) and (6). As to initial condition, u has the above 
tanh profile, v = 0, T = 1, p = 1, p = (7 Af 2 ). -1 . The 
Reynolds and Prandtl numbers are 174000 and 0.707, re- 
spectively. The nozzle radius R, which is the reference 
length /*, is taken to be 1.95 cm. 

A 600 x 160 stretched grid, which extends up to 
about 59.3 and 9.2 radii in the respective axial and ra- 
dial directions, is used. Figure 2, drawn at every fourth 
grid line, shows the grid clustered near the nozzle height. 
Instantaneous contour plots of the vorticity, Mach num- 
ber, and static pressure are presented in Figure 3. These 
clearly indicate a large vortical structure in the devel- 
oping jet flow. Local pressure minima coincide with the 
centers of vortices. Figure 4 is the time sequence of the 
vorticity of the jet. Each frame of the contour plots is 
100 A* apart. The flow structure is wavy yet well con- 
nected in the early stage of the flow development region. 
This wave motion is magnified and the confined vorticity 
layer rolls up to form a discrete vortex lump, which can 
be seen in every frame of vorticity contour. This vortex 


roll-up is followed by vortex shedding and vortex pairing 
as the flow proceeds downstream. The vortex pairing is 
the coalescence process of two consecutive vortices and 
is known as the shear layer growth mechanism. The con- 
vective velocity with which the large vortical structure 
moves is obtained to be 1.06 from Figure 4, which is 
very close to the mean velociy of the two streams. The 
convective velocity is constant regardless of the loca- 
tion and size of the vortex. The only exception is when 
the vortex pairing takes place where the vortex which 
travels behind the preceeding one speeds up to catch up 
with the latter and the two slide on each other to coalese 
into a larger one. The convective velocity obtained by 
the present computation agrees very well with the ex- 
pressions given by Dimotakis 14 and Papamoschou and 
Roshko 15 . They assumed that the dynamic pressure of 
the two streams should be about the same in the frame 
which moves with the convective velocity. 

Figure 5 shows the axial velocity spectra. The 
abscissa is Stf(D), the Strouhal number defined to be 
fD/Uj with U j the jet center velocity at the nozzle 
exit and the ordinate represents the absolute value of 
Fourier coefficient. D is the diameter of the nozzle. The 
most preferred frequency of St(D)— 0.63, which appears 
dominant at the upstream region, persist up to x=4 D 
in Figure 5. This dominant high frequency shifts to- 
wards the low frequency and the spectra show the single 
most dominant frequency of about St(D)= 0.35 from the 
downstream of x=6jD. This frequency is close to the av- 
erage Strouhal number of 0.3 based on the downstream 
puff counts by Crow and Champagne 16 . Many experi- 
ments have confirmed that the St(D) of about 0.3 is the 
dominant frequency of the organized vortical structure 
where the potential core ceases to exist. 

Peaks in the power spectra of the u velocity show 
the most preferred frequencies : St(#):=0.006 at x—Dj 2 
and jD, St(0)=O.O21 at x=2D, St(0)=O.O22 at x=3 D, 
and Stf(0)=O.O34 at x=4 D. This Strouhal number is de- 
fined to be f0/Uj with 6 the local momentum thickness 
defined as : 

0= A \p - U >)Wf- U ) dr (8) 

St(6) is the most frequently used flow parameter in de- 
scribing initial flow development. According to the lin- 
ear theory by Michalke 17 , the most amplified St(0) is 
about 0.017 for a spatially evolving jet flow with a hy- 
perbolic tangent mean velocity profile. The linear the- 
ory agrees with Freymuth’s experiment 3 up to about 
5t(^)=0.01, but experiment exhibits a flat peak area 
over 0.01 < St(9) < 0.025. Considering a flow region in 
which flow solution adjust itself from a boundary condi- 
tion to a flow solution, the computed Strouhal number 
of about 0.02, which is most preferred at x=2D and 3D, 
is close to both the theory and the experiment. 
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The mean axial velocity profiles are given in Figure 
6. The potential core in the middle erodes as flow pro- 
ceeds downstream and vanishes as shown at x=8 D. Fur- 
ther downstream the mean velocity profile maintains the 
same profile as at r=8 D rather than undergoing rapid 
decay. This is a huge departure from the experimental 
observation. Figure 7 also shows the profiles of u rm s- 
The growth of tirms the nozzle height (i.e. r=l) is 
shown in Figure 8. Both natural and forced cases are 
presented. Forcing is given at the most preferred up- 
stream St(D) of 0.63 and its subharmonic 0.315. u rms 
becomes saturated at about 0.3 for all three cases. The 
maximum of u r ms grows up to about 0.3 and never de- 
cays far downstream. An Euler computation shows es- 
sentially the same results as the viscous computation. 
This is because our viscous computations do not include 
an adequate device dissipating turbulence energy due to 
small scale random motion. 

Far-field Sound Prediction 

As mentioned earlier, the flow field as simulated will 
be used as the sound source. There are two main ways to 
predict far-field noise from the numerical flow solution. 
The first approach is the acoustic analogy. According to 
LighthilFs acoustic analogy, the wave equation for static 
pressure p can be obtained as : 


d 2 p 

d& 



d 2 

dxidxj 


(pU.Uj - Tij) 


( 9 ) 


where a 0 is the reference ambient speed of sound. A 
deviation from an isentropic state, which is d 2 {p — 
a lp)/dt 2 i is neglected in equation (9). The source 
term in the right-hand-side of the above equation is of 
quadrupole type. If the source term is assumed a pri- 
ori known and physical obstacles are not present, the 
solution of equation (9) is given by a volume integral as 

p(x,<)-P 0 = ^ j v \[qY dV(y) (10) 

where r = |x — y|, po is the undisturbed pressure, q the 
source terms in equation (9) and the quantity in [ ]* is 
defined as : 

[ 9 ]*=9(y,< + ) = g (y,t-^) 


t+ is called the retarded time, which is the emission time 
for a acoustic signal to travel from the source point y 
to the field point x. Figure 9-(a) shows the geometric 
variables used in the above definition. P(x) is the field 
point where p is being sought and Q(y) is the source 
point. More details of the acoustic analogy method can 
be found in Fuchs and Michalke 18 . Direct computa- 
tion of the far-field sound by the acoustic analogy from 


the unsteady flow solution can be found in the refer- 
ences[19,20]. The acoustic analogy is a useful method 
to predict a far-field pressure. However, because it uses 
the wave equation as a governing equation even in the 
flow region where all kinds of nonlinear interactions oc- 
cur, the solution omits important wave phenomena such 
as the refraction due to shear and the doppler effect by 
movement of the source element. 

If the wave phenomenon is separable from the fluid 
dynamics, the wave equation cam be solved to take into 
account all the physical effects. Let us imagine that 
away from the flow area there is a region where flow 
fluctuation is negligible so that a pure acoustic field as- 
sumption is valid. In this case, the flow and the acoustic 
regions can be separated and the flow equation and the 
wave equation can be solved side by side in their respec- 
tive regions. Acoustic information is transfered from the 
flow solution on the boundary joining the two regions. 
This surface away from a flow regime is illustrated in 
Figure 9-(b). If all the necessary pressure information is 
provided on the surface, the pressure at x can be given 
by the Kirchhoff surface integral form as : 


p(x,0 -Po 


47T 



1 dr 
do r dn 


[|n«(y) 


( 11 ) 


where [ ]* has the same definition a is before and n the 
unit normal to the surface. Equation (11) is the solution 
of the wave equation for the acoustic medium at rest. 
The necessary information includes the pressure and its 
normal and time derivatives on the surface. 

Wave Solution Example 

A pressure wave, which is generated by a simple 
monopole source and radiated spherically, has been cho- 
sen as an example problem to demonstrate the surface 
integral method given by equation (11). The pressure 
satisfying a homogeneous form of the wave equation (9) 
can be written as : 


p — p 0 = sin 2n(r - t)/r (12) 

Length and time are nondimensionalized by wave length 
X and inverse of frequency / _1 . We now choose a cylin- 
der surface as an arbitrary boundary surface on which 
the histories of p,p n ,Pt are provided. This is shown in 
Figure 10. The source is located at the origin. The 
cylinder radius is 0.5 and length is unity, so the cylinder 
is not acoustically compact. For a spatial resolution in 
the surface integral, Az and A r are taken to be 3^. Cir- 
cumferntial elements are set to be 100 on the side and 
two lid surfaces of the cylinder. For retarded time con- 
sideration, a complete set of pressure information over 
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an entire time period is needed. Far-field points, which 
are observation points, are located 10° apart at a dis- 
tance r = 50A as shown in Figure 10. Only angles from 
0° to 90° are considered due to symmetry. Figure 11 
plots the sound pressure (i.e. pnns) against observa- 
tion angle and shows how crucial it is to resolve the 
retarded time. Exact value of the prms for this example 
is 0.01414 with no preferred directivity, prms obtained 
by employing 40 time elements is almost the same as 
the exact solution. Figure 12-(a) also shows the time 
history of the pressure at angles 0° and 90°. The two 
curves show little difference when 40 time segments are 
used in resolving a time period. It can be concluded 
that a proper resolution of the retarded time is essential 
to get the correct wave characteristics. Figure 12-(b) 
is for an instantaneous pressure against distance which 
decays a s r - ” 1 . Two curves at 0° and 90° are identical 
with 40 time elements per period. 

If the source has multiple frequency contents, the 
refinement of the retarded time is increasingly impor- 
tant since the highest frequency of interest should be 
adequately resolved. Likewise, the spatial refinement 
also becomes important. Therefore, for wave propaga- 
tion, which is generated by a distributed source with 
frequencies varying over a broad range, refinement of 
surface elements should be emphasised as well. 

Far-Field Sound Generated by Unsteady Jet Flow 

Figure 13 gives the schematics of the surface on 
which the pressure information is specified and the field 
positions where the noise is observed. A cylinder is cho- 
sen to accomodate a jet column. The radius of the cylin- 
der is about 4 R and the length is about 20 R. The obser- 
vation points are on the plane comprising the axis of jet 
center spacing 10° apart circumferentially at distances of 
200 R and 400it!. For the frequency of about S*(D)=0.3, 
which is fluid dynamically dominant downstream, an 
acoustic wave length can be estimated to be about 11 
for the jet of M=0.6 speed. This length could be much 
shorter for the wave with upstream preferred Strouhal 
number. Therefore, the cylinder body chosen can not 
be treated as acoustically compact, so the retarded time 
contribution should be considered. The length of time 
record we have to keep track of depends on the size of 
the cylinder. The difference between the minimum and 
maximum times to transmit an acoustic signal to the 
observation point is about 16 (i.e. times taken for the 
sound to travel from A to B and from A to C). With 
At of 0.00875, 1867 time steps are needed to cover the 
retarded time which spreads over 16 time units. This is 
too long, so the acoustic computation is performed at 
every 3 time steps. We set A t s — 3A t and 630A t s are 
used to retarded time distribution in the noise compu- 
tation. The At and A t s are referred to as the flow and 


acoustic time steps, respectively. 

The cylinder surface is divided into 80 circumferen- 
tial elements, but only half of them are used due to the 
axisymmetry. In performing the surface intergal of equa- 
tion (11), a contribution by two lid surfaces is omitted 
and only the side surface of the cylinder is considered. 
258 mesh points in the axial direction are used to form 
a cylinder which extends to x=20i? downstream with 
the radius of 4 R. Therefore, an additional storage of 
258 x 630 for each of p,p n , and p t and 258 x 40 x N 0 
for each of r and dr/dn is needed. N 0 is the number of 
observation points, which is 19 in the present work, and 
40 is the half of the circumferential elements. 

The OASPL (the overall sound pressure level) at 
r=200i? and 400 R is given in Figure 14. The com- 
puted result qualitatively agrees with the experiment by 
Moore 21 which revealed the maximum in the OASPL 
at about 20° over a variety of subsonic speed condi- 
tions. The differences are that the present OASPL has 
the maximum at about 10° and double-peak in 10° to 
40° , whereas the experimental data shows a single peak 
in this area. Figure 15 shows the spectra of the pressure 
at r=100D. At 0°, peak Strouhal numbers are 0.13, 
0.20, 0.28, 0.46. At 10°, peak Strouhal numbers are 
0.13, 0.20, and 0.46. From the pressure spectra between 
0° and 20° the dominant frequencies fall in the range 
of 5<(jD)=0.13~ 0.28. As the angle diverges from the 
flow direction, the dominant peak shifts toward higher 
frequencies, which is Strouhal number of 0.61 briefly at 
30° and of 0.46 for a broad range of angle. 

Conclusion and Discussion 

Far-field sound pressure has been evaluated directly 
from a side-by-side computation with the solution of un- 
steady jet flow. The OASPL has maximum directivity 
at about 10° from the jet flow axis. The spectra of 
the far-field pressure show a similar trend as the exper- 
iment in that low frequency dominates at the smaller 
angle from the flow axis and shifts to the higher fre- 
quency as the angle increases. Although the present re- 
sults agree qualitatively with the experimental observa- 
tions, the present method has the following limitations 
: (1) helical mode of the jet flow is excluded due to 
the axisymmetry assumption, which has a strong effect 
on the growth of the instability downstream. (2) The 
turbulence intensity never decays even far downstream 
because no adequate dissipating mechanism is provided 
for small scale turbulence energy. (3) The Green’s fun- 
tion used in the Kirchhoff surface integral for the wave 
is best suited to wave propagation problems through a 
quiescent medium, so the surface integral is to be modi- 
fied to be the solution of the convective wave equation if 
the surrounding fluid has some considerable flow speed. 

Items (1) and (2) suggest that the flow simulation 
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be of LES type so that the jet flow solution can accomo- 
date a non-axisymmetric mode and mimic a realistic de- 
cay of the jet column downstream. This would probably 
lead to a quantitatively more accurate noise prediction. 
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Figure 1. Wave relation given by Equation (4). 
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Figure 3. Contour plots of (a) vorticity, (b) Mach number, and (c) static pressure. 
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Figure 5. Spectra of the u velocity. 


Figure 7. urms profile. 
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Figure 15. Spectra of the pressure at the distance of 100D. 
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